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ABSTRACT 

Context, Procyon is one of the brightest stars in the sky and one of our nearest neighbours. It is therefore an ideal benchmark object 
^3 ■ for stellar astrophysics studies using interferometric, spectroscopic, and asteroseismic techniques. 

r-j ' Aims. We used a new realistic three-dimensional (3D) radiative-hydrodynamical (RHD) model atmosphere of Procyon generated 

with the Stagger Code and synthetic spectra computed with the radiative transfer code Optim3D to re-analyze interferometric and 
I . spectroscopic data from the optical to the infrared. We provide synthetic interferometric observables that can be validated against 

O ' observations. 

Methods. We compute intensity maps from a RHD simulation in two optical filters centered at 500 and 800 nm (Mark III) and one 
infrared filter centered at 2.2 fim. (Vinci). We constructed stellar disk images accounting for the center-to-limb variations and used 
them to derive visibility amplitudes and closure phases. We computed also the spatially and temporally averaged synthetic spectrum 
from the ultraviolet to the infrared. We compare these observables to Procyon data. 

Results. We study the impact of the granulation pattern on center-to-limb intensity profiles and provide limb-darkening coefficients 
in the optical as well as in the infrared. We show how the convective related surface structures impact the visibility curves and 
closure phases with clear deviations from circular symmetry from the 3rd lobe on. These deviations are detectable with current 
interferometers using closure phases. We derive new angular diameters at different wavelengths with two independent methods based 
on 3D simulations. We find 6y m d = 5.390 ± 0.03 mas and prove that this is confirmed by an independent asteroseismic estimation 
(^seismic = 5.360 ± 0.07 mas. The resulting T sff is 6591 K (or 6556 K, depending on the bolometric flux used), which is consistent with 
Teff.iR = 6621 K found with the infrared flux method. We find also a value of the surface gravity logg = 4.01 ± 0.03 [cm/s 2 ] that is 
larger by 0.05 dex from literature values. Spectrophotometric comparisons with observations provide very good agreement with the 
spectral energy distribution and photometric colors, allowing us to conclude that the thermal gradient of the simulation matches fairly 
£Nj , well Procyon. 

Finally, we show that the granulation pattern of a planet hosting Procyon-like star has a non-negligible impact on the detection of hot 
Jupiters in the infrared using interferometry closure phases. It is then crucial to have a comprehensive knowledge of the host star to 
directly detect and characterize hot Jupiters. In this respect, RHD simulations are very important to reach this aim. 
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1. Introduction ISchaeberlddl896l). It was one of the first detected white dwarfs 



(Eg gen & Greensteinll965l) . The main component of the system 
Procyon (a Cams Majons) is one of the brightest stars in the sky is a subgiant F5 IV - V (Procyon A, HR 2943, HP 61421) which 



and one of our nearest neighbours. It is therefore an ideal target ise ndin g its li fe on the Main Sequence (lEggenberger et al.l2005 

for stellar astrophysics studies. [ Provost et all l2006h. It has a sol ar metaliicitv dOriffinl 1 1 97 1 , 

For this reason, it has a long history of observations. |Bessel| |Steffenlll985b Ullende Prietol^ 12001 with an e ffective tem- 



<ll844i) discovered that its motion was perturbed by a invisible pe rature around T„« * 6500 K rtCode et aLlf!?^ . One of the 



companion. Procyon became, after Sinus, one of the first astro- first radius determinations was made photometrically by iGravl 

metric binar ies ever detecte d. The first orbital elements were de- d!%7h . who found R = 2.24 R Q . The binary nature of the system 

termined by|Auwers| dl862|), who showed that the period of rev- is a great opportunity to determine t he mas s of bo th compan- 

olution is about 40 years. The faint companion, Procyon B, was ions xhe first attempt was made by Istnmdl (fl95lh who deter- 
not detected visually until the end of the nineteenth century by 
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mined the masses of Procyon A & B, by determining the par- 
allax and orbital elements of th e system. He f ound 1.76 ±0.1 
and 0.65 + 0.05 M Q , respectivelv. lSteffenl d 19851) claimed that the 
mass of Procyon A was too large to be compatible with its lu- 
minosity and suggested rather a mass around 1.4 M Q . More re- 
cent and more accurate determinations of the mas s of Procyon 
A con verge to a star with M = 1.497 + 0.037 M n dGirard et all 
120001) or 1.430 + 0.0034 M Q dGatewood & Hanll2006l) . The age 
of the system is well constrained by the white dwarf compan- 
ion, whose cooling law as function of time is well established. 
iProvencal et alj d2002l) found an old white dwarf with an age of 
1.7 + 0.1 Gy. This determination is a strong constraint for stellar 
evolution models. 

A way to discriminate between different masses and ages 
is the determination of the interferometric radius. The first 
attempt to measure t he diameter of Procyon was done by 
Hanbury Bro wn et all (fT%7l fl97l who found an angular di- 
ameter of 6 — 5.50 + 0.17 mas . This value was confirmed 
later by iMozurkewich et aD d 19911) but with a much b etter pre- 
cision (~ 1%). More recently. iKervella et al] d2004bl) redeter- 
mined the angular diameter using the Vinci instrument at VLTI. 
They found a diameter that is even smaller 9 = 5.448 + 0.053 
mas. lAufdenberg et alj d2005l) re-analyzed these data using hy- 
drodynamical model atmospheres and found 6 = 5.404 ± 0.031 
mas (see their Table 7). It is interesting to compare these angular 
diameters with the independent i nfrared flux method, a r ecent 
determination was performed by ICasagrande et al.l (l2010h who 
derived a value of (fa = 5.326 ± 0.068 mas, which is smaller 
than the lKervella et all s result, but they agree to within la. 

Another great particularity of Procyon is the presence of os- 
cillations that are due to trapped acoustic modes. The first claims 
of a detectio n of an excess power w ere made by iGellv et all 
dl986lll988l) and lBrownet ail (119911) who found a mean large 
spacing between consecu tive acoustic mode s of about 39 //Hz 
and 55 //Hz, respectively. IMartic et al.l (1 19991) made a clear de- 
tection using the ELODIE fiber echelle spectrograph with a 
strong excess power aroun d v max » 1000 //Hz and confirmed 
the results of iBrown et al.l dl99ll) . They determined the fre- 
que ncy spacing to be 5 5 //Hz. Later, lEggenberger et alj d2004l) 
and IMartic et alj (|2004) made the first identifications of indi- 
vidual frequencies of spherical harmonic d egrees I = 0,1,2 
with mean large spacing s of 55.5 ± 0-5 (Eggenbergeret al.l 
120041) and 53.5 ± 0.5 //Hz dMartic et al.ll2004l) . More recent ob- 



servations have been made from the ground during single or 



multisite campaign s dMosser et al ] 12008b lArentoft et al.l 12008: 
IBedding et al.lboiOl) or from space on board of t he MOST satel 



lite dMatthews et al.|[2004T: iGuenther et al.ll2008l) . which consid- 
erably improved the precision in frequencies within ~ 1 //Hz. 

The stellar evolu tion model of Procyon was made by 
iHartmann et al.1 d 19751) . who used the astrometric mass, photom- 
etry, and the angular diameter of lHanburv Brown et al.1 (fl97l 
to constrain the model. More realisti c modeling, especially 
with b etter equation-of-states, came with lGuenther & Demarqua 
d 1993b . who showed the importance of diffusion of elements. 
The values of individual p-modes frequencies considerably 



(Barban et al. 


19991 Idi Mauro & Christensen-Dalsgaardl 1200 ll 


Eggenberger et al. 


120051: IProvost etaTI 120061 iBonanno et al.l 



20071: IGuenther etal.1 120081) . lEggenberger et al 



and 



Provost et al. d2006l) made a realistic stellar evolution modeling 



to constrain simultaneously the lo cation in the HR d i agram and 
the large frequency separations of lEggenberger e t al. (2004J) and 
IMartic et al.l ( 20041). respe ctively. They disagreed on the derived 
mass. Eggenberge r et al.1 found a mass that agrees well with 



the astrometr ic value derived by iGirard et al.l d2000l) . whereas 
IProvost et al.l found a mass that c orresponds to the more recent 
value of iGatewood & Hani d2006l) . Part of the difference in the 
two stellar evolution models comes from the slightly different 
large separation (which increases with the mass), and partly due 
to the difference in the stellar evolution codes themsel ves. We 
emphasize that the best stellar model of IProvost et al.l that fits 
asteroseismic and spectro-photometric data has an age (> 2 Gy) 
that is consistent with the age of the white dwarf. Indeed, the 
astrometric mass of Procyon B implies that the mass of its pro- 
genitor was about 3M Q , which has lived about ~ 500 My on 
the Main Sequence. Therefore, we can safely conclude that the 
age of the system must be at least 2 Gy. The larger mass of 
lEggenberger et al.l d2005l) corresponds to a younger (< 1.7 G y) 
system, incompatible with the result of IProvencal et al.l d2002l) . 

The atmospheric parameters (T & ff,logg, [Fe/H]) and the in- 
terferometric radius, that are used to define the stellar evolu- 
tion model and the analysis of frequencies, strongly depend 
on the realism of the atmosphere and the exactness of the 
temperature gradient in the surface layers. In the case of F- 
stars, these gradients are strongly modified by the convec- 
tive transport which is more vig orous than in the Sun. Thi s 
is clearly seen in line bisectors dG7a^fl98l |D ravin s1 fT987h . 
which are two or three times larger than in the Sun. The 
higher convective velocities are due to the higher stellar lu- 
minosity and smaller densities. This strong effect of convec- 
tion must be taken into account for stellar physics diagnos- 
tics. Realistic 3D time-dependent hydrodynamical simulations 
of the s urface layer of Procyon were performed by several 
authors dAtroshchenko etanil989t iNordlund & Dravinsl [19901 
lAllende Prieto et al.l 12002 ). who showed that the 3D effects in 
such F-star can bring significant differences for line profil e 
formation and abu ndance analysis ( Allende Prieto et al.ll2002l) . 
iNelsonld 19801) and lNordlund & Dravinsl d 19901) showed that the 
surface defined at T = T e g is not flat but rather "corrugated" 
due to the large fluctuations and the high contrast of granula- 
tion. These hydrodynamical simulations can reproduce with suc- 
cess the line shifts, asymm etries, and in particular the observed 
bisect ors of various lines dAllende Prieto et al 1 120021 Fel and 
FeII). lAilende" Prieto et ail also showed that the 3D limb darken- 
ing law significantly differ ID law, up to ~ 1.6%, leading to a 
correction of Ar e ff « 50 K, which is non-neg l igible for precise 
stellar evolution modeling. lAufden b erg et al.l d2005l) made 3D 
models using the CQ 5 BOLD code dFrevtag et alj|2002ll201 lb to 
calculate limb darkened intensity profile s to analyze visibility 
curves obtained by the Vinci instrument dKervella et alj|2003bl) 
and the Mark III. T heir 3D analysis led to a smaller radius by 
0.04 mas tha n that oflKervella et al.1 d2003bl) . who used ID limb 
darkened law. lAufdenberg et al.1 also showed that 3D models bet- 
ter reproduce the spectral energy distribution in the UV, whereas 
ID model are unable. 

Regarding the importance of a realistic hydrodynamical 
modeling of the atmosphere of Procyon, we re-analyze the in- 
te rferometric and spectrosc opic data at the different wavelengths 
of lAufdenberg et all d2005l) using up-to-date line and continuum 
opacities dGustafsson et al.l2008l) to derive a new radius. We pro- 
pose a solution that agrees well with the astrometric, asteroseis- 
mologic, infrared flux method, and interferometric data. We also 
explore the impact of convection-related surface structures on 
the closure phases and assess how the direct search of planets 
by interferometry may be affected by the host star surface struc- 
tures. 
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2. Three-dimensional radiative-hydrodynamical 
approach 

2.1. Procyon simulation 

The c onvective surface of Pr ocyon is modeled usingthe Stagger 
Code dNordlund et al.ll2009l Nordlund & GalsgaarcQ). In a local 
box located around the optical surface r » 1, the code solves 
the full set of hydrodynamical equations for the conservation of 
mass, momentum, and energy coupled to an accurate treatment 
of the radiative transfer. The equations are solved on a staggered 
mesh where the thermodynamical scalar variables (density, inter- 
nal energy, and temperature) are cell centered, while the fluxes 
are defined on the cell faces. This scheme has several numerical 
advantages to simulate surface convection. It is robust against 
shocks and ensures conservation of the thermodynamic vari- 
ables. The domain of simulation contains the entropy minimum 
located at the surface and is extended deep enough to have a 
flat entropy profile at the bottom. The code uses periodic bound- 
ary conditions horizontally and open boundaries vertically. At 
the bottom of the simulation, the inflows have constant entropy 
and pressure. The outflows are not constrained and are free to 
pass through the boundary. The code is based on a sixth order 
explicit finite difference scheme and fifth order for interpolation. 
Numerical viscosity of the Rytchmeyer & Morton type is used to 
stabilize the code. The corresponding adjustable parameters are 
chosen to minimize the viscosity and are not adjusted to fit the 
observables. We used a realistic equation-of-state that accounts 
for io nization, recombination, and dissociation dMihalas et al.l 
1 1988b and continuous (T rampedach et al. privat e communica- 
tion) and line opacities (Gustaf sson et al.1 120081) . An accurate 
treatment of the transfer is needed to get a correct temperature 
gradient of the transition region between optically thin and thick 
layers. The transfer equation is solved using a Feautrier-like 
scheme along several inclined rays (one vertical, eight inclined) 
through each grid point. The wavelength dependence of the ra- 
diative transfer is taken into account using opacity bins lNordlundl 
(Il982h . The numerical resolution is 240 3 . The geometrical sizes 
are 22 Mm x 22 Mm horizontally and 17 Mm vertically. The 
horizontal dimensions of the box are defined to include at each 
time step a sufficient number of granules and the vertical one is 
chosen to ensure that the entropy profile is flat at the bottom. The 
equations of magnetohydrodynamics are not computed for this 
model. The stellar parameters corresponding to our RHD model 
(TableQ} are r eff = 65 1 2 ±25 K, log g = 4.0 [cm/s 2 ], and a solar 
chemical composition (lAsplund et al. 120091) . The uncertainty in 
T e ff represents the fluctuations with time around the mean value. 
These parameters roughly correspond to those of Procyon. The 
exa ct values of the parameters do not influence the limb darken- 
ing (lAufdenberg et al.ll2005l:lBigot et al.ll201lt) . 

2.2. Spherical tiling models, intensity maps, and spectra 

The computational domain of each simulation represents only 
a small portion of the stellar surface. To obtain an image of 
the whole stellar disk, we emplo y the same tiling method ex- 
plained in lChiavassa et al ] (l2010al) . For this purp ose, we used the 
3D pu re-LTE radiative transfer code Optim3D (IChiavassa et alj 
120091) to compute intensity maps from the snapshots of the 
RHD simulation of Table[T]for different inclinations with respect 
to the vertical, fi= cos(6»)=[ 1.000, 0.989, 0.978, 0.946, 0.913, 
0.861, 0.809, 0.739, 0.669, 0.584, 0.500, 0.404, 0.309, 0.206, 
0.104] (Fig. [TJ and for a representative series of simulation's 

1 1995, http://www.astro.ku.dk/~kg/Papers/MHDcode.ps.gz 



M = 0.40 




u h i-0 ii 



Fig. 1. Intensity maps in the Mark 500 nm filter for the RHD 
simulation for different inclination angles /i. The intensity ranges 
from lxlO 3 to 1.1 x 10 7 erg cm 2 s _1 A -1 . 



snapshots: we chose ~ 25 snapshots taken at regular intervals 
and covering ~1 h of stellar time, which corresponds to ~5 p- 
modes. Optim3D takes into account the Doppler shifts due to 
convective motions and the radiative transfer equation is solved 
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monochromatically using pre-tabulated extinction coefficients as 
functions of temperature, density, and wavelength (with a resolv- 
ing power of A/6A = 500 000). The lookup tables were computed 
for the same chemic al composition as the RHD simulation (i.e. 
lAsplund et al.ll2009b and using the same extensive atomic and 
mol ecular opacity data as t he latest generation of MARCS mod- 
els dGustafsson et al]l2008l) . 

Then, we used the synthetic images to map onto spherical 
surfaces accounting for distortions especially at high latitudes 
and longitudes cropping the square-shaped intensity maps when 
defining the spherical tiles. Moreover, we selected intensity 
maps computed from random snapshots in the simulation time- 
series: this process avoided the assumption of periodic boundary 
conditions resulting in a tiled spherical surface globally display- 
ing an artifactual periodic granulation pattern. 

Based on the stellar radius estimates and on the sizes of the 
simulations domains (Table [TJ, we required 215 tiles to cover 
half a circumference from side to side on the sphere (number of 
tiles = it ■ R+/22.0, where 22.0 is the horizontal dimension of the 
numerical box in Mm and the radius of the star). To produce 
the final stellar disk images, we performed an orthographic pro- 
jection of the tiled spheres on a plane perpendicular to the line- 
of-sight (// = 1.0). The orthographic projection returned images 
of the globe in which distortions are greatest to ward the rim of 
the hem isphere where distances are compressed dChiavassa et al] 
l2010ah . 

In this work, we computed a synt hetic stellar disk image fo r 
interferometric s pectral bands use d in lAufdenbere et al.l (I2005h : 
(i) the Mark I II dShao et al.ll 19881) c entered at 500 and 800 nm, 
and (ii) Vinci dKervella et al.ll2003al) centered at 2.2 /mi (Fig. [2]). 
The Mark III sensitivity curves are assumed to be Gaussian with 
central wavelengths Ap = 500 and 800 nm, ea ch of them with 
a FWHM of 20 nm dMozurkewich et al.lfl99"ll) . We produced a 
number of synthetic stellar disk images corresponding to differ- 
ent wavelengths in the filters with a spectral resolving power of 
20 000. Figure[3]shows the resulting synthetic stellar disk images 
averaged over each passband. 

With Optim3D, we computed also the spectra, which are nor- 
malized to the filter transmission as: ^^7777-77-, where I A is the 

J T{A)dA 

intensity and T (A) is the transmission curve of the filter at a cer- 
tain wavelength. The spectra in Fig.|2]have been computed along 
rays of four //-angles [0.88, 0.65, 0.55, 0.34] and four 0-angles 
[0°, 90°, 180°, 270°], after which we performed a disk integra- 
tion and a temporal average over all selected snapshots. 



2.3. Three-dimensional limb-darkening 

Figure [3] shows irregular stellar surfaces with numerous convec- 
tive related surface structures. There are pronounced center-to- 
limb variations in the Mark 500 nm and Mark 800 nm filters 
while these are less noticeable in the Vinci filter. This is manly 
due to the different Planck functions in the optical range and in 
the infrared region. 

We derived azimuthally averaged (i.e., averaged over differ- 
ent <p angles) intensity profiles for every synthetic stellar disk 
i mage from the simulation (Fig.H] ). Using the method described 
in Chiavassa et al. I d2009l l2010alh the profiles were constructed 
using rings regularly spaced in // = cos(fl) for fj. < 1 (i.e. 
r/R* < 1), with 9 the angle between the line of sight and the 
vertical direction. The standard deviation of the average inten- 
sity, (Tkjj), was computed within each ring, the // parameter is 
connected to the impact parameter r/R* through the relation- 
ship r/R* = -y/l - 11 2 , where R* is the stellar radius reported in 
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Fig. 2. A spatially and temporally averaged synthetic spectrum 
in the optical wavelength region with the Mark 500 nm filter 
(red: top left), Mark 800 nm (green: top right), and Vinci in- 
frared filter (blue: bottom). 



Table Q] The total number of rings is 20, we ensured that this 
number is enough to have a good characterization of the inten- 
sity profiles. 

Figure |4] displays a steeper center-to-limb variation for the op- 
tical region, as already visible in the disk images, with fluctua- 
tions of ~ 20% in the Mark III 500 nm filter down to ~ 10% 
and ~ 5% in the Mark III 800 nm and Vinci filters, respectively. 
We tested the impact of spectral lines in the Mark 500 nm filter, 
for which the effects of lines are stronger, computing a synthetic 
disk image considering only the continuum opacities. Its aver- 
aged intensity profiles (red line in left panel of Fig. |4| is very 
similar to the one computed with spectral lines (black line) with 
differences lower than 1% for r/R* < 0.9 and between 1 - 5% 
at the limb (0.9 < r/R* < 1.0). It is also interesting to notice 
that the continuum profile tends to be closer to the uniform disk 
profile (dashed line in the figure), as well as the Mark 800 nm 
and Vinci filters with respect to the intensity profile of the Mark 
500 nm. 

W e used the following limb-darkening law dChiavassa et al.l 
120091) to fit the averaged profiles: 



(i) 



k=0 
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Table 1. 3D simulation of Procyon used in this work. 



< r eff >l [Fe/H] 


logg 


x, y, z-dimensions 


x, y, z-resolution 


R* 


[K] 


[cm/s 2 ] 


[Mm] 


[grid points] 


[Ro] 


6512+25" 0.0 C 


4.0 


22.0x22.0x17.0 


240x240x240 


2.055" 



" Horizontally and temporal average and standard deviation of the emergent effective temperatures 

b ICollet etal.N201lh 

c Chemical composition by Asplund et al. (2009) 

d Angular diameter of 5.443 mas (Kervella et al. 2004b]) converted into linear radius with Eq. ifTJ 



Mark 500 nm Mark 800 nm Vinci 




, [Rsun] x [Rsun] » [Rsun] 

Fig. 3. Top row: synthetic stellar disk images of the RHD simulation. The intensity range is [0.3 - l.OxlO 6 ], [0.3 - 5.2xl0 5 ], and 
[0.3 - 0.9xl0 5 ] ergcrrT 2 s _1 A -1 for Mark 500 nm, Mark 800 nm, and Vinci filters, respectively. Bottom row: enlargements of the 
images above. 

Mark 500 nm Mark 800 nm Vinci 




Fig. 4. Radially averaged intensity profiles (black line) derived from the synthetic stellar disk images of Fig. [3] The gray areas 
denote the one sigma spatial fluctuations with respect to the averaged intensity profile. The dashed horizontal line is an uniform 
disk intensity profile. The intensity is normalized to the mean intensity at disk center and the radius is normalized to the radius 
given in Table Q] The blue dotted-dashed line in central and right panels is the average intensity profile of Mark 500 nm. The red 
dotted-dashed line in left panel is the continuum only averaged intensity profile in Mark 500 nm. 
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Mark 500 nm 




Table 2. Limb-darkening coefficients from the law described in 
text for RHD simulation of Table [T] 



0.0 



Mark 800 nm 



4 




0.0 



Vinci 



4 




0.0 



Fig. 5. Center-to-limb fits (red dashed line) obtained using 
Eq. ([TJ with N — 2 for the 3D azimuthally average intensity 
profile (solid line). A full limb-darkening (dash-dotted line), a 
partial limb-darkening (triple dot-dashed line) are also shown, fi 
= cos 6, with 9 the angle between the line of sight and the vertical 
direction. 



Filter 


a 


a i 


fl2 


Mark 500 nm 


1.000 


-0.066 


-0.421 


Mark 800 nm 


1.000 


-0.041 


-0.236 


Vinci 


1.000 


0.016 


-0.188 



solution with very minor improvement s to the y minimization 
using N — 3. This was already found bv lChiavassa et al.l(l2010al) 
for K giants. Figure [5] shows the fits in the different filters used 
in this work, while Table [2] reports the limb-darkening coeffi- 
cients from the fits. The figure illustrates that the average inten- 
sity profiles from a RHD simulation is largely different from the 
full limb-darkening (/(//)//( 1) = jj.) and partial limb-darkening 
(/(//)//( 1) = 0.5 • yu) profiles as well as the power law profile that 
does not even provide an appreciable fit to the radially average 
profile. We therefore discourage the use of these simple laws. 



2.4. Visibility curves 

The synthetic disk images are used to derive interferometric ob- 
servables. For this pu rpose, we used the method described in 
IChiavassa et al.l d2009t) to calculate the discrete complex Fourier 
transform T for each image. The visibility, vis, is defined as 
the modulus of the Fourier transform normalized by the 
value of the modulus at the origin of the frequency plane, \To\, 
with the phase tantp = 3(T)/1R(T), where 3(T) and %(T) 
are the imaginary and real parts of the complex number f , re- 
spectively. In relatively broad filters, such as Vinci, several spa- 
tial frequencies are simultaneously observed by t he interferom- 
eter. Th is effect is called bandwidth smearing. iKervella et al.l 
(l2003bllcl) show that this effect is negligible for squared visibil- 
ities larger than 40% but it is important for spatial frequencies 
close to the first minimum of the visibility function. To account 
fo r this effect, we compute d the squared visibilities as proposed 
bv lWittkowski et all (12004 



f vis A A 
{vis ) = — , 



(2) 



where vis\ is the squared visibility at wavelength A, Aq and 
A\ are the filter wavelength limits, is the transmission curve 
of the filter, and Fx is the flux at wavelength A (see Fig. [2j». vis\ 
has been computed from the disk images that already include 
in the calculations. 
In our case we follow these steps: 

1 . we generate a synthetic stellar disk image at different wave- 
lengths of the filters; 

2. we compute the visibility curves vis\ for 36 different cuts 
through the centers of the stellar disk images; 

3. we apply Eq. (O to obtain the averaged squared visibilities. 



In this equation, h(p) is the intensity, a* are the limb- 
darkening coefficients, and N + 1 their number. We performed 
a Levenberg-Marquardt least-square minimization to fit all the 
radially averaged profiles of Fig. fusing this law and weighting 
the fit by l/cr /(K) , due to 3D fluctuations (gray areas of Fig. |4j 
We varied the order N and found that N - 2 provides the optimal 



Figure [6] (top panels) shows the visibility curves computed 
with the Eq. (0 and for 36 different cuts through the centers 
of synthetic disk images. This is equivalent to generating dif- 
ferent realizations of the stellar disk images with intensity maps 
computed for different sets of randomly selected snapshots. A 
theoretical spatial frequency scale expressed in units of inverse 
solar radii (Rq 1 ) is used. The conversion between visibilities ex- 
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pressed in the latter scale and in the more usual "arcseconds" 
scale is given by 

vis ["] = vis [Rq 1 ] • d [pc] • 214.9 (3) 

where 214.9 is the astronomical unit expressed in solar ra- 
dius, and d is the distance of the observed star. The spatial fre- 
quency in arcsec -1 (i.e, v) is related to the baseline (i.e., E) in 
meters by 



where A is the wavelength in /urn. 
The first null point of the visibility is mostly sensitive to 
the ra dial extension of the observ ed object (e.g. lOuirrenbachl 
1200 U and IChiavassa et all 1201 Obi for an application to RHD 
simulations), while the first null point and the second lobe 
of the visibility curves are sensitive to the limb-darkening 
dHanburv Brown et al.lll974h . Since we want to concentrate on 
the small scale structure of the surface, the visibility curves of 
Fig. [6] are plotted longward of the first null point. They show 
increasing fluctuations with spatial frequencies due to devia- 
tions from the circular symmetry relative to uniform disk visi- 
bility. This dispersion is clearly stronger in the optical filter at 
500 nm and appreciable from longward of the top of third lobe. 
Moreover, it is also noticeable that the synthetic visibilities are 
systematically lower than the uniform disk with a weaker di- 
vergence for the Vinci filter. This is due to a the non-negligible 
center-to-limb effect visible in Figs. [3] and [4] 

The bottom panel of Fig.|6]displays the one cr visibility fluc- 
tuations, F, with respect to the average value vis (F = cr/vis). 
The dispersion increases with spatial frequency. This is due 
to the small scale stru cture on the model stellar disk (see e.g. 
IChiavassa et al 1201 Obi) . The dispersion is stronger in the case of 
the Mark III 500 nm filter with respect to the redder filters. 
Figure|4]shows small differences between the averaged intensity 
profiles of the Mark III 500 nm filter with and without consid- 
ering the spectral lines in our calculations. Therefore, we com- 
puted visibility curves in both cases and found that the visibility 
fluctuations are indistinguishable (Fig. |6l bottom panel). While 
the molecular absorption can cause a strong difference in stel- 
lar surface appearance (and consequentially also on the visibil- 
ity curves) in the case of cool evolved stars (e.g., the contribu- 
tion of H2O to the radiu s measurement for red supergiants stars, 
Chiavassa et al 1 l2010bh . this is not the case for Procyon where 
the atomic lines are not strong enough to cause an appreciable 
effect and also the atmosphere is very compact. 
The bottom panel of Fig. [6]also shows that, on the top of the sec- 
ond lobe (~ 0.4 Rq 1 ), the fluctuations are of the order of ~ 0.5% 
of the average value for Mark III 500 nm filter and ~ 0.1% for 
Mark III 800 nm and Vinci filters. From the top of the third lobe 
(~ 0.6 Rq 1 ) on, the fluctuations are ~ 2% of the average value 
in the optical region, whic h is larger of the instru mental error of 
VEGA on CHARA (-1%, iMourard et alll2009h . 
It must be noted that our method of constructing realizations 
of stellar disk images inevitably introduces some discontinu- 
ities between neighbouring tiles by randomly selecting temporal 
snapshots and by cr opping intensity maps at high latitudes and 
longitudes. However. lChiavassa et al.l(l2010al) proved that for the 
signal artificially introduced in to the visibility curves is largely 
weaker than the signal due to the inhomogeneities of the stellar 
surface. 



3. Fundamental parameters 

3. 1 . Multiwavelength angular diameter fits in the Mark /// 
and Vinci filters 

We used th e data from the Mark I II interferometer at 500 
and 800 nm (iMozurkewich et alJll99lh plus additional measure- 
ments at 50 nm reported in lAufdenberg et all d2005l) . It has to 
be noted that lMozurkewich et alj d!99 lb do not provide a calibra- 
tor list even thought they claim that the calibrators are all smaller 
that the science star. The Mark III 500 nm filter has points on 
the first and the second lobe while the Mark III 800 nm has lim- 
ited baselines in the first lobe. For the Mark III 500 nm filter, we 
used only the data with baselines larger than 20 meters because 
at shorter baselines the error bars are > 10% of the squared visi- 
bility. We also used Vinci data from t he Very Large Telescop e 
Interferometer dKervella et alj l2004at lAufdenberg et al.ll2005l) . 
The wavelength of these observations is around 2.2 ^m and the 
baselines go up to 66 m, significantly lower in the first lobe of 
the Procyon visibility curve. 

We used two independent methods to determine the angular 
diameters: 

1 . the synthetic visibilities of Fig. |6]obtained from the spherical 
tiling method described in this work and scaled to absolute 
model dimensions (using Eq. (O and (|4]i) to match the inter- 
ferometric observations; 

2. the method described in lAufdenberg et al.l d2005l) : 
iBigot et all d2006l 1201 lb . that computes the normalized 
fringe visibilities using the van Cittert-Zernike theorem. 

Method (1) is a the direct consequence of the Fourier 
Transform of the synthetic disk images (Fig. [3): the visibil- 
ities vary for different cuts through the centers of synthetic 
disk images (i.e., the position of the first null changes its po- 
sition). Method (2) is based on the integration of the spatial av- 
erage intensity profile. For both methods, we used a Levenberg- 
Marquardt least-square minimization. 

Table[3jreports the best-fit angular diameters. While for the uni- 
form disk there is a significantly wide range of values among 
the different filters, the angular diameters from the RHD sim- 
ulation are closer and overlap with the uncertainties. The two 
independent methods used in Table [3] show a clear tendency 
of the optical diameters to appear smaller than in the infrared 
( fl g" k "' ~ 0.99). The possible explanation for the different ra- 
dius between Mark III and Vinci are: (i) the RHD model used is 
not fully appropriate to the observations, for instance the center- 
to-limb variation may not be correct because, for instance, the 
points on the second lobe (sensitive to the center-to-limb varia- 
tion) of Mark 500 nm filter do not match very well (Fig. [7] top 
panel); (ii) the data of Mark III may present problems (e.g., cal- 
ibrator or systematics). 

In addition to this, we also conclude that the two methods give 
consistent results and can be used without distinction to perform 
angular diameter fits at least for dwarf stars. If fact, in the case 
of cool evolved stars with low surface gravity (logg < 2.0), the 
surface asymmetry may strongly impact the shape of the star 
and the radius depend s on orientation of the projected baseline 
dChiavassa et al]l2009h . 

The angular diameters found for the Vinci filter (Tab le [3]) 
are in fairly good agreement with lAufdenberg et al.l (|2005) who 
found 5.403 + 0.006 mas. However, gy inci is ~ 2% smaller than: 

(i) 5.448 + 0.03 mas dKervella et al.ll200 4b). based on the fit of 
the Vinci with baseli ne points of Fig. [7] lower than 22 meters; 

(ii) 5.48 ± 0.05 mas dAllende Prieto et al]l2002l) . obtained with 
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Fig. 6. Top panels: visibility curves (vis) from the Procyon simulation in the filters computed in this work. The visibilities are 
computed for 36 different azimuth angles 5° apart (thin gray lines). The dotted line is an uniform disk scaled to match the same 
radius, i.e. the first null in visibility curve. The visibilities are displayed only longward of the first null visibility point. A logarithmic 
scale is used on y-axis. Bottom panel: visibility fluctuations (a /vis) with respect to the visibility average value (vis) as a function of 
spatial frequencies for all the considered filters. 



Table 3. Best angular diameters, in mas, for RHD simulation 
with Method (1) (£>Methodi) an d Method (2) (£>Method2) described 
in Section |3~T1 and for the uniform disk. The error is one <x with 
respect to the average value. 



Filter 


$Methodl 


#Method2 


VD 




Mark 500 nm 


5.313 ± 


0.03 


5.324 


± 0.03 


5.012 ± 


0.03 


Mark 800 nm 


5.375 ± 


0.06 


5.370 


+ 0.06 


5.208 ± 


0.06 


Vinci 


5.382° i 


0.03 


5.397" 


±0.03 


5.326 ± 


0.03 



" Since the resulting diameters are very similar, we adopted an aver- 
age angular diameter of 5.390 for all the calculations in the article 



Eq. [7] fro m the linear radius of 2.07 1 ± 0.020 R Q ; (iii) 5.50 ± 
0.17 mas dMozurkewich et all 199l|). Our Ov^-, is also close to 
5.326 ± 068 mas found bv ICasagrande et all d2010h with the 
infrared flux method. 

Transforming the Vinci angular diameter to linear radius 
(Eq.[7]), we found 2.019 R . This value is compared to what we 
chose as a reference in TablefT] 2.055 R . The difference is 0.023 
R Q , which is 16 Mm (i.e., 0.7 times the numerical box of RHD 
simulation, 22 mM, Table [T] 4th column). We checked that this 
has a negligible impact on the observables found in this work 

using the spherical tiling method described in Section | 2.2l 

The astrometric mass (M = 1.430 + 0.034 M jGatewood & Hani 
[20061 used in Section 13.3b combined with our interferomet- 
ric diameter leads to a new gravity \ogg = 4.01 + 0.03 
[c m/s 2 l, which is larg e r by 0.05 dex than the value derived 
in lAllende Prieto et all ( 120021) . The contribution of the revised 
Hipparcos parallax is ~0.01 dex in the log g change. 



3.2. The effective temperature 

The effective temperature is determined from the bolometric flux 
Fbol and the angular diameter 9 in the Vinci filter (average be- 
tween the two methods in Table[3] 9 = 5.390 mas) by the relation 



T e f[ - 



AF 



BOL 



CT st 9 2 



0.25 



(5) 



where <x st stands for the Stefan-Boltzma nn constant. Our de- 
rived d ia meter being smaller by ~ 2% than iMozurkewich et al.l 
dl991l) or lKervella et alj d2004bl) . the deriv ed effective tempera- 
ture ha s to be larger by ~ 1 % compa red to lAllende Prieto et al] 
(120021) or lAufdenberg et all (120051) . who used the previous 
mentioned references. There are several sources for the bolo- 
metric flux leading to slightly differe nt T e g. The values of 
Fbol = (18.20 + 0.43) 10~ 6 erg cirT 2 s _1 dFuhrmann et al.ll 19971) 
and F bql = (17.82 + 0.89) 10~ 6 erg cm' 2 s" 1 dAufdenberg et al I 
120051) lead to r eff , 3D = 6591 ±43 K and 6556±84 K, respectively. 

Our new 3D 7 W returns a value clos er to T e g ^ = 6621 ± 
80 K obtained by ICasagrande et alj d2010l and Casagrande pri- 
vate communication) tha n the old derived value of 6516 ± 87 K 
(lAufdenberg et al.ll2005l) . 

The influence of the uncertainties in the selected fundamental 
parameters (r e jf,logg,[Fe/H]) of our RHD model atmosphere 
has a negligible impact on the limb-darkening and therefore 
on the deri ved angular diam eter and T e ff. This was tested for 
HD 49933 dBigot et al.ll201lh . a star similar to Procyon. 
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0.447 




Fig. 7. Best matching synthetic visibility from Fig. [6] for the 
Mark 500, Mark 800, and Vinci data. The stellar disk maps in 
each filter have been scaled to match the radius of the star. For 
the Mark 500 nm filter, we used only the data with baselines 
larger than 20 meters. 



3.3. Asteroseismic independent determination of the radius 

The radius of the star can be derived from its oscillation spec- 
trum. The frequency v max of the maximum in the power spectrum 
is generally ass umed to scale with the acoustic cut-off frequenc y 
of the star (e.g. iBrown etal 1 I199U iKieldsen & Beddingll201 lh . 



Table 4. Photometric colors for RHP simula t ion of Table [TJ and 
for the corresponding observations of lBessell (1 19901) 





B- V 


V-R 


R-I 


V-I 


RHD simulation 


0.419 


0.255 


0.252 


0.507 


Observation 


0.420 


0.245 


0.245 


0.490 



therefore v„ 
the radius 



R 

~rZ 



gl VTeff- Then, it is straightforward to derive 



M v n 



Mr, 



0.5 



' eff.O 



-0.25 



(6) 



The validity of such scaling rel ation has been verif i ed on 
large asteroseismi c surveys (e.g. iBedding & Kieldsenl 12003b 
IVerner et al.l 1201 lT) . The value of v max derived from photom 
etry i s accurately determined, 1014 + 10 //Hz ( Arentoft et al 
120081). The solar value of v max is taken from Belkace m et al 
(1201 ll) . We emphasize that the dependence on T e g in Eq. [6] is 
weak, therefore the derived radius is not very sensitive to the 
selected value of the effective temperature. We use the value of 
T e ff = 6591 ± 43K derived in Section I3~2l since it is closer to 
the infrared flux method determination. Since Procyon is a bi- 
nary star, the mass can by determined by the astrometric orbital 
elements and the third Kepl er's law. However, t he derived value 
is subject to debate. Indeed. iGirard et aTl ( 20001) found a m ass of 

found 



Gatewood & Ha 



M = 1 .497 + 0.037 M Q , whereas|< 
M = 1 .430 + 0.034M o . As discussed in the introduction, we pre- 
fer to keep the value of iGatewood & Hani since stellar evolution 
models that use this mass, are in agreement with the age of the 
white dwarf companion. Using these values, we found a radius 
of R = 2.023 + 0.026 Rq. We can translate this radius into angu- 
lar diameter using the relation 



Seismic = 2 (R/R Q ) n p tan(e Q /2) = 5.36 + 0.07 mas, 



(7) 



with the solar ang ular radius 9 Q /2 = 959.64 + 0.02 arcsecs 
(IChollet & Sinceacjl 19991) and the parallax n p = 284.56 ± 1.26 
mas (Ivan Leeuwenl2007l) . This radius agrees well with our inter- 
ferometric value within error bars. 



4. Spectrophotometry and Photometry 

The radiative transfer code Optim3D includes all the up-to-date 
molecular and atomic line opacities. This allows for very real- 
istic spectral synthetic RHD simulations for wavelengths from 
the ultraviolet to the far infrared. It is then possible to compute 
realistic synthetic colors and compare them to the observations. 
Figure [8] shows the spectral energy distribution (SED) computed 
along rays of four /i-angles [0.88, 0.65, 0.55, 0.34] and four (p- 
angles [0°, 90°, 180°, 270°], and after a temporal average over 
all selected snapshots (like in Fig. |2j. The shape of SED reflects 
the mean thermal gradient of the simulations. 

We used the prescriptions by iBesselll (1 19901) to compute the 
color indexes for the filters BVRI. Table |4]displays the compar- 
ison of the synthetic colors with the observed ones. The differ- 
ence is very small. 

We used the absolute spectrophotometry measure- 
ments at ultravi ol et an d visual wavelengths collected by 
lAufdenberg et al.l (120051) to compare with the synthetic 
SED. The data come from: (i) the Goddard High Resolution 
Spectrograph (GHRS) data sets Z2VS0105P ( PI A. Boesgaard), 
Z17X020CT, Z17X020AT, Z17X0208T (PI J. Linsky) from 136 
and 160 nm; (ii) the International Ultraviolet Explorer (IUE) 
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Fig. 8. Spatially and temporally average synthetic spectral en- 
ergy distribut ion from the U V to the near infrared. Red curves 
correspond to lBesselll d 19901) for the filters BVRI. 



iRodnguez-Pascual et al ] (119991) from 170 to 306 nm; (iii) the 
Hubble Space Telescope imaging Spectrograph (STIS) from 220 
to 41 nm; (iv) and the visua l and near-infrared wavelength data 
from lGlushnevaet all (Il992l) . The data from the GHRS are far 
superior to any other measurements below 160 nm because the 
continuum drops by more than a factor of 100 here, too much 
for the limited dynamic range of IUE. The flux was estimated 
by computing the mean flux between the emission lines in each 
spectrum incorporating the flux uncertainties provided with 
each calibrated data set ( lAufdenber g et al. 20051). 

Figure [9] (top row) shows that our synthetic SED matches 
fairly well the observations from the UV to the near-IR region. 
This result can be compared to the already good agreement 
fo und with composi te RHP model done with CO 5 BOLD code 
in lAufdenberg et al.l (120051) . Moreover, the bottom row of the 
Figure shows that at higher spectral resolution, the agreement 
is even more remarkable. RHD models match also the obser- 
vations between 136 and 160 nm, that are supposed to for m at 
depths beneath its chromosphere (lAufdenberg et al]l2005l) and 
thus being affected by the convective-related surface structures. 
It should also be noted that the ultraviolet SED longward of 
160 nm may be impacted by non-LT E treatment of iron-group 
elements. IShort & Hauschildl d2005l) found that non-LTE mod- 
els for the Sun have up to 20 % more near-UV flux relative to 
LTE models. It is not possible to determine if this difference is 
also present in Procyon because either the RHD model calcu- 
lation and the post-processing calculations have been done with 
LTE approximation. The ultraviolet SED may also be affected by 
scattering at these wavelengt hs. This effec t is als o not included 
in our calculations. However, iHaveketaU (120101) demonstrated 
that, in RHD simulations, the scattering does not have a signifi- 
cant impact on the photospheric temperature structure in the line 
forming region for a main sequence star. 

We conclude that the mean thermal gradient of the simulation, 
reflected by the spectral energy distribution, is in very good 
agreement with Procyon. 



5. Closure phases and perspectives for hot Jupiter 
detection 

Interferometry has the potential for direct detection and char- 
acterization of extrasolar planets. It has been claimed that dif- 
ferential interferometry could be used to obtain spectroscopic 
information, planetary mas s, and orbit inclination of extrasolar 



planet s around nearby stars: ISegransan et al.l d2000l); Lopez et al.l 
(l2000h: iJoergens & Ouirrenbachl (120041) : iRenard et all (|2008|) : 
Zhao et alj d2008l) : Ivan Belld d2008l) : iMatter et alf J2OIOI) : 
Absil et all d201lh : IZhao et al.l d201 ll) . However, current inter- 



ferometers lack sufficient accuracy for such a detection. When 
observing a star with a faint companion, their fringe patterns 
add up incoherently and the presence of a planet causes a slight 
change in the phases and, consequently, the closure phases. This 
difference can be measured with a temporal survey and should 
be corrected with the intrinsic closure phases of the host star. 
In this framework, theoretical predictions of the closure phases 
of the host stars are crucial. Closure phase between three (or 
more) telescopes is the sum of all phase differences: this proce- 
dure removes the atmospheric contribution, leaving the phase in- 
formation of the object visibility unaltered. The major biases or 
systematic errors of closure phases come from non-closed trian- 
gles introduced in the measurement process, which in principle, 
can be precisely calibrated. Therefore, it is a good obser vable 
for stable and precise measurements (e.g. lMonnier|l2007l) . The 
closure phase thus offers an important complementary piece of 
information, revealing asymmetries and inhomogeneities of stel- 
lar disk images. 

5. 1 . Closure phases from the star alone 

Figure [10]displays deviations from the axisymmetric case (zero 
or ±ji) that are particularly occurring in optical filters. There is 
a correlation between Fig. [10] and Fig. |6]because the scatter of 
closure phases increases with spatial frequencies as for the vis- 
ibilities: smaller structures need large baselines to be resolved. 
Moreover, it is visible that the closure phase signal becomes im- 
portant for frequencies larger than ~ 0.6 R ' (top of the third vis- 
ibility lobe, i.e., vis ~ 0.06 from Fig. |6j. These predictions can 
be constrained by the level of asymmetry and inhomogeneity of 
stellar disks by accumulating observations on closure phase at 
short and long baselines. 

Observing dwarf stars at high spatial resolution is thus cru- 
cial to characterize the granulation pattern using closure phases. 
This requires observations at high spatial frequencies (from 3rd 
lobe on) and especially in the optical range. In fact, consid- 
ering that the maximum ba seline of CHARA is 331 meters 
dten Brummelaar et alj | 2005[). for a star like Procyon with a par- 
allax of 284.56 mas dvan Leeuwenl l200~7h . a baseline of ~55 
(240) meters is necessary to probe the third lobe at 0.5 (2.2) /urn. 
For comparison, the nominal measurement error with CHARA 
array is 0.3° with a peak of performance of 0.1° for a shorter 
triangle IZhao et al] d2008l) (to be compared with closure phases 
value of Fig. [Toll. 

5.2. Closure phases from the hosting star plus the hot 
Jupiter companion 

Among the detected exoplanets, the direct detection and the 
characterization of their atmospheres appears currently within 
reach for very close planet-star system (<0.1AU) and for plan- 
ets with temperatures >1000 K, implying their infrared flux is 
~ 10~ 3 of their host stars. Since the bulk of the energy from 
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Fig. 9. Comparisons between synthetic SED (black) and spectrophotometric measurements at ultraviolet and visual wavelengths 
(green dots with error bars in left column plots and red line). Top row: the models and data are bin ned to 2 nm resolution i n the UV 
for clarity and to 5 nm in the visual to match the resolution of the observed spectrophotometry of iGlushneva et al] dl992l) . Bottom 
row: same as above with the synthetic SED and the observations at higher spectral resolution. The synthetic SEDs are scaled in 
absolute flux using an angular diameter of 5.390 mas found from the fit of Vinci filter observations. 



hot Jupiters emerges from the near-infrared between 1-3 /jm 
(iBurrows et al.1 120081) . interferometry in the near infrared band 
(like in the range of the Vinci filter centered at ~ 2.2 yum) can pro- 
vide measurements capable to detect and characterize the plan- 
ets. Detecting hot Jupiters from Earth is endeavoring because it 
is challenging and at the limits of current performance of inter- 
ferometry and, once these conditions are met, the signal from 
the host stars must also be known in detail. As visible in the 
synthetic stellar disk images of Fig. [3] the granulation is a non- 
negligible aspect of the surface of dwarf stars and has an impor- 
tant signal in the closure phases (Fig.fTOii. 

To estimate the impact of the granulation noise on the planet 
detection, we used the RHD simulation of Procyon and added a 
virtual companion to the star. The modelling of the flux of an ir- 
radiated planet requires careful attention on the radiative transfer 
conditions rel ated to the stel l ar irra diation and therefore we used 
the models of iBarman et al.l (1200 ll) . In particular, we used spec- 
tra of hot irradiated extrasolar planet around a star with about the 
same spectral type of Procyon, a mass of 1 Jupiter mass, and an 
intrinsic temperature 1000 K. We assumed a radius of 1 .2 Jupiter 
radii a nd various orbital d istances [0.1, 0.25, 0.5, 1.0] AU fol- 
lowing iMatter et al.l(l2Q10h . The atmospheric composition of the 



hot Jupiter is identical to the two models of lAllard et al.1 (1200 ll) 
where: (1) dust (particles and grains) remains in the upper at- 
mosphere and (2) where dust has been removed from the upper 
atmosphere by condensation and gravitational settling. 

Figure QT| shows the geometrical configuration of the star- 
planet system for a particular distance. First we average the ex- 
oplanetary spectrum in the range of the Mark III 800 nm and 
Vinci filters (we had no data for the exoplanet spectrum in the 
range of Mark III 500 nm filter), and then we used this inten- 
sity for the stellar companion in the Fig.[TT|. The intensity of the 
planet is stronger in the infrared with respect to the optical. The 
ratio between the stellar intensity at its center (i.e., p = 1), /vinci 
or /Mark800> an d the planet integrated intensity in the same filter, 

^planet? IS. 



'planet 



= 11.5,31.7,42.7,46.3, 



jMarkgOO 
^planet 



34.7, 104.0, 130.0, 182.5 



(8) 
(9) 



for distances [0.1, 0.25, 0.5, 1.0] AU, respectively. 



11 



Chiavassa et al.: Three-dimensional interferometric, spectrometric, and planetary views of Procyon 



Mark 500 nm 



Vinci 



200 



100 7 



■100 



-200 



200 



100 



■ ■ . »■ ■»'.y*iW » *i l K ' < 



0.5 



1.0 



1.5 



2.0 



Ro" 1 




Mark 800 nm 



-100 7 



-200 



200 



100 



Fig. 11. Synthetic stellar disk images (yellow-orange scale) of 
the simulation in Vinci filter (the intensity range is [0.3 - 
0.9xl0 5 ] erg cirT 2 s _1 A -1 ) together with a hot Jupiter (pink) at 
a distance of 0.1 AU for a star with a parallax of 284.56 mas 
(Ivan Leeuwenl l20O7h . The intensity of the hot Jupiter is ~ 1 1 
times lower than the maximum stellar surface intensity. 
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Fig. 10. Scatter plot of closure phases of 2000 random baseline 
triangles with a maximum linear extension of ~ 155 m at 500 
nm, ~250 m at 800 nm, and ~680 m at 2.2 n m, respectively. 
A parallax of 284.56 mas (Ivan Leeuwenll2007f) and an apparent 
radius of 5.390 mas have been used. 



We computed the closure phases from the image of Fig.lTTI 
and for similar systems corresponding to the Vinci and Mark III 
800 nm filters and for the distances reported above. These clo- 
sure phases were compared with the resulting phases computed 
for exactly the same triangles but for a system without the pres- 
ence of a planet. Figure [12] displays the absolute differences be- 
tween closure phases with and without the presence of a hot 



Jupiter. Setting as a reference the closure phase nominal error 
of CHARA (0.3°) and also 1° (horizontal lines in the plot), only 
the Vinci filter gives differences that should be detectable on the 
third lobe (0.5 < R e l < 0.8) while there is no signature on the 
second lobe (R e l < 0.5). The absolute difference increases as a 
1 5 20 function of spatial frequency. It is indistinguishable for all the 

planet' distances for frequencies larger than 0.8~ R^ 1 . 

The purpose of this Section was to show the alteration of 
the signal due to the granulation noise in the detection of a hot 
Jupiter around Procyon-like stars using closure phases. So far, 
no companions have been detected around Procyon except for a 
white dwarf astrometric companion detected already in the 19th 
century (see Section[TJ. 

The studies of hot Jupiter atmospheres will reveal their com- 
position, structure, dynamics, and planet formation processes. It 
is then very important to have a complete knowledge of the host 
star to reach these aims. 



6. Conclusions 

We have provided new predictions of interferometric and 
spectroscopic observables for Procyon, based on RHD hydro- 
dynamical simulation, that affect the fundamental parameter 
determination of the star and are important for the detection hot 
Jupiter exoplanets. 

We have studied the impact of the granulation pattern on the 
center-to-limb intensity profiles and provided limb-darkening 
coefficients in the optical as well as in the infrared. We showed 
that synthetic visibility curves from RHD simulation are sys- 
tematically lower than uniform disk and this effect is stronger 
in the optical filters. In addition to this, visibilities display 
fluctuations increasing with spatial frequencies (i.e., departure 
from circular symmetry) that becomes ~ 2% on the top of third 
lobe in Mark III 500 nm filter. However, in the Mark III 800 
nm and Vinci filter the dispersion is much weaker. 

We have derived new angular diameters at different wave- 
lengths with two independent methods based on the RHD 
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Fig. 12. Absolute differences of the closure phases computed for 
the 2000 random baseline triangles in the case of a system with- 
out the presence of a hot Jupiter companion (CP/2 nop i an et) and a 
system star-planet (CPh p u net , see text and Fig. [Till. The horizon- 
tal solid lines indicate the 1.0 and 0.3 degrees and the vertical 
dashed line the spatial frequency corresponding to the end of 
third visibility lobe. The black, blue, green, and red colors de- 
note the planet-star distance of [0.1, 0.25, 0.5, 1.0] AU, respec- 
tively. The maximum linear extension, the apparent radius and 
the parallax are the same as Fig.lTOl 

simulation. The angular diameter in the Vinci filter is 9 = 5.390 
mas, which leads to an effective temperature of T e ff = 6591 ± 43 
K or 6556 ± 84 K, depending on the bolometric flux considered. 
This value is now more consiste nt with r e ff.rR - 6621 ± 80 
K from the infrared flux method dCasagrande et al.l I2010L and 
Casagrande private communication). 

Using an independent estimation of the radius from asteroseis- 
mology, we found se i S mic = 5.36 + 0.07 mas. This radius agrees 
well with our interferometric value within the error bars. 
Eventually, the combination of the astrometric mass and 
our new interferometric diameter leads to a new gravity 
logg = 4.01 + .03 fcm/s 2 ], which is larger by 0.05 dex than the 
value derived in lAllende Prieto et al.l (l2002h . 

We provided accurate comparison of synthetic spectrum 
from a RHD simulation to observations from the ultraviolet to 
the infrared. The photometric BVRI colors are very similar to 
the observations while it is difficult to conclude on the infrared 



colors because of the saturation of the observations. Also the 
comparison with t he absolute spec t ropho tometric measure- 
ments, collected in I Aufdenberg et alj d2005l) . is in agreement. 
We conclude that the mean thermal gradient of the simulation, 
reflected by the spectral energy distribution, is in very good 
agreement with Procyon. 

The convective related surface structures impact also the sig- 
nal of the closure phases that show departures from symmetry at 
about the same spatial frequencies of visibility curves. We con- 
cluded that closure phases not equal to or +n may be detected 
with today interferometers such as CHARA in the visible filters 
where the baselines are long enough to get to the second/third 
lobes. 

We estimated the impact of the granulation noise on the hot 
Jupiter detection using closure phases around stars with the same 
spectral type of Procyon. We used the synthetic stellar disks ob- 
tained from RHD in the infrared and optical filters and added 
a virtual companion to the star based on real integrated spectra 
of irradiated extrasolar planet. Then, we computed the closure 
phases for planet-star system and star only and found that there 
is a non-negligible and detectable contamination to the signal of 
the hot Jupiter due to the granulation from spatial frequencies 
longward of the third lobe. This is valid only for the infrared 
where the energy from the hot Jupiters is stronger. It is thus very 
important to have a comprehensive knowledge of the host star 
to detect and characterize hot Jupiters, and RHD simulations are 
very important to reach this aim. In a forthcoming paper, we will 
extend this analysis to solar type stars and K giants across the HR 
diagram. 
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